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A systematic evaluation turn boon conducted o£ 
certain digital imago onlmncomant techniques per- 
formed in Image uimce. Three types of images hnvo 
boon uucd, computer gcnorntud plirnttoms, tomograms 
of n ayntbocic plnmtom, mu! axial tomograms of 
human anatomy containing Imagca of Icaionu, arti- 
ficially introduced into tlio tomograms . 

Several typoa of smoothing, sharpening, and 
histogram modification have been explored. It has 
been concluded that cite moat useful enhancement 
techniques aco n selective smoothing of singular 
picture elements, combined with contrast manipula- 
tion, Tho moot useful tool in applying these tech- 
niques is tho grayscale histogram. 

Stimmn rv 

NASA developed technology lias heen applied to 
the enhancement of computerised axial tomograms of 
human anatomy, Tho work was performed on an image 
analysis system developed for tho Inspection of 
acroopnce structures. Tho tomograms wore supplied 
by tho Clovelsnd Clinic Foundation in tho form of 
cumputar readable magnetic tapes. 

A systematic evaluation was conducted of « 
repertory of digital enhancement techniques per- 
formed in image space. Three typos of images wore 
used! computer gonorntod phantoms, tomograms of a 
synthetic phantom, and tomograms of human anatomy 
containing images of lesions, artificially Intro- 
duced Into the tomogramB. 

The repertory of enhancement functions Included 
uuveral types each of smoothing, sharpening and 
histogram manipulation. The functions were applied 
in an orderly manner to several of each of the 
three typos of images. The input parameters were 
varied systematically, 

The following conclusions were readied for the 
typo of imagos investigated in tills programs 

1, Selective smoothing of singular picture ele- 
ments is most useful for reducing seise spikes. 

2, Contrast manipulation is the most effective en- 
hancement method, 

3, The grayscale histogram is the most useful 
• tool in applying these techniques , 


The author extends his gratitude to Thomas F, 
Heaney, M.D., Chairman, Division of llndlology, and 
Jeffry M. Arnold, t'h.D. , (formerly) Head, scienti- 
fic Computing, of tho Cleveland Clinic Foundation, 
for providing tho tomograms uboJ in this work. A 
special word of thanks is duo to Mr. Sanford F, 
Felder, Technology Utilisation Engineer, NASA Lewis 
Research Center, for coordinating the program. 

Introduction 

Tito raw images produced by computerised axial 
tomography possess certain shortcomings. Fine 
structure is frequently obscured by limited reso- 
lution and noise introduced by both the X-ray sen- 
sors and the reconstruction algorithm. These short- 
comings are shared witli thane of radiogrsphic images 
produced in the nondestructive inspection of ma- 
terials and structures. The resolution of such prob- 
lems in aerospace applications la currently under 
investigation at NASA, Lewis Research Center. A 
real-time computerized radiographic inspection fa- 
cility lias been developed^, 

To explore common problems, a cooperative pro- 
gram between tho Nondestructive Evaluation Section 
of Lewis Research Center and tho Division of Radi- 
ology of the Cleveland Clinic Foundation lias been 
initiated, It is coordinated by tho Technology 
Utilization Office of Lewis Research Center, Tho 
work is being conducted by NASA personnel in tho 
laboratories of .ho Nondestructive Evaluation Sec- 
tion. 

Axial tomograms of human anatomy arc supplied 
by Cleveland Clinic. They arc transferred in the 
form of computer readable magnetic tape. They are 
processed on Uiu NDE section computer, using a 
reportory of enhancement functions dovolopcd for 
on-line radiographic inspection. Currently, tho 
investigation is limited to operations in image 
space. It is planned to be extended to Fourier 
filtering mid nnulyeis in Fourier space. 

The tomographic images supplied by Cleveland 
Clinic consist of 256 lines of 256 square picture 
elements. The imago analysis hardware processes 
images of 490 linos of 512 rectangular picture ele- 
ments, Tho elements luivo nti aspect ratio of 4 j 3, 
tho broadcast television standard. For display 
purposes, it is necessary to transform the Cleveland 
Clinic images by cyclically repeating picture ele- 
ments. However, they are stored in tho computer 
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times memory, and processed, as 256 X 256 arrays. 
The imago analysis hardware ubos a range of 
brightness values from 0 (black) to 255 (white). 
The tapes supplied by Cleveland Clinic use arbi- 
trary X-ray absorption values, They ave trans- 
formed on reading to "delta numbers" from -1000 
to represent the X-ray film density of air (black) 
through +1000 to represent bone (white). 

Procedure 


Three classes of enhancement functions were 
explored: image smoothing, image sharpening, and 

histogram modification. These techniques ore des- 
cribed more fully elsewhere*. Four types of image 
smoothing functions were used: 

1. Weighted gray value averaging over a neigh- 
borhood, 3 picture elements wide by 3 high. The 
weight given each element is on input parameter. 

2, A weighted median otherwise identical to item 

1 . 


3, A replacement of a neighborhood of picture ele- 
ments by the mean of the neighborhood. The size 
may be made 2X2, 4X4, 3X8, or 16 X 16. 

4, Selective averaging, the criterion being that 
the central point in a 3 X 3 neighborhood be out 

of a specified range. If the element is not out of 
range, it is not averaged. 

Two types of sharpening functions were explored: 

1, Lap la clan, 

2. Root-moan-square gradient. 

Two types of histogram modification were inves- 
tigated: 

1. Contrast expansion, in which a specified input 
range of gray values is transformed to a different 
specified output range. Picture elements falling 
outside the input range are transformed to the 
minimum or maximum value. 

2. Gray value segmentation in which an arbitrary 
number of gray value boundaries is specified, and 
any given input value is transformed to the mean 
of the range in which it falls, 

These functions were systematically applied to 
3 types of images: 

1. Computer generated phantoms having irregular 
boundaries. The boundaries, gray values, and 
superimposed noise are input parameters. A typi- 
cal phantom image had an average gray value on one 
side of the boundary of 100, and 104 on the other 
side (in terms of the hardware values). Random 
noise having a range of ±2 was superimposed on the 
image, producing a range of 98 to 102 on one side, 
□nd 102 to 106 on the other. This image, contrast 
expanded to full scale, is shown in figure 1. 



Figure 1, - Computar generated phantom. 

2. Several tomograms of a synthetic phantom. 

This was a nylon body of irregular shape immersed 
in varying strength solutions of potassium iodide 
and water. One of these images is shown, contrast 
expanded to full scale, in figure 2. 



Figure 2. Tomogram of synthetic phantom. 


3. Several tomograms of human unatomy, to which 
images of lesions were artificially introduced . 

One of these is shown full scale in figure 3, with 
the locations of the artificial lesions indicated 
by the labels a, b, c, d, e, f, and g. 

In all cases, the anomolies were not Immediate- 
ly apparent to the unaided human eye. 

Several sets of input parameters were tested 
far each function. Each function was, in turn, 
applied to several of the 3 types of images. Com- 
binations of functions were used when that appeared 
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Figure ). - Tomogram uf Imnun anatomy. Pointers 
Indicate hu nt Ion of U'n Ioiih . 


fo he a logical procedure. The output Images were 
Inspected visually In an attempt to locate the 
known bounder leu or ananollea. When neceaaary, 
histograms of the vartoua images were produced. 

The gray-value hiatograra, preaented on the video 
monitor, la a plot of the number of picture ete- 
menta corrcapondlng to each gray value. It la taken 
over a rectangular region defined by a curaor dls- 
played on the monitor. Typically, It covera ap- 
proximately 40 percent of the total Image area. It 
la ahown aa a white rectangle In figure 3. The 
hlatogram la normalized to a maximum ordinate value 
of 100. The number of picture elements at 0 (black) 
may be omitted from the normalization computation 
In response to a orompi by the computer program. 

The plot Illustrated In figure 4 uses the brlght- 


I* . 1 t 



O 


IOO 


I'LK'tLMI > u 

nFnwh - i?o 

RHTO iV°- 

NOX COUNI Ol 


i'OO 

uo 


Kl U' inrtOL. I ULL UCflLL 


ness values produced by ths Image analysis hard- 
ware (0 to 2)1). 

The histograms were used both to define tha 
values of Input parameters to the various functions, 
and to analyze their effect on the output Image. 

Investigations conducted on the computer gen- 
erated phantoms revealed that random noise reduces 
the effectiveness of the two types of Image sharp- 
ening functions. In the absence of noise, both of 
these functions clearly reveal the boundary be- 
tween areas of differing gray value, even when they 
are virtually Indistinguishable by the unaided eye. 

A thresholding process Is especially effective. 

In this case, a histogram Is generated of the 
gradient Image. The histogram Is used to define a 
threshold value. All values greater than this 
threshold are set to black, and all others set to 
white. This binary Image Is then logically summed 
with the original Image, often with dramatic re- 
sults. 

However, the presence of random noise intro- 
duces a speckle pattern superimposed over the gra- 
dient Image. This pattern can be minimized either 
bv smoothing the original image, or by selectively 
discarding isolated black picture elements on the 
binary gradient Image. But both procedures tend to 
obscure the boundary. When the range of noise Is 
such that there Is an overlap of gray values be- 
tween the two regions, the methods fail completely. 
This is the case with the computer generated phan- 
tom described herein, and la apparently the case 
with the tomograms. 

The various smoothing procedures reduce the 
Isolated noise spikes In the inwiges. However, the 
more effective the procedure In reducing noise, the 
more It obscures boundaries. The only one of these 
functions which leaves recognizable boundaries Is a 
selective averaging procedure. In this procedure, 
each picture element Is compared with the range of 
Its 1X3 neighbors. If It falls outside a de- 
fined multiple of that range, It Is replaced by the 
mean of the neighborhood. Otherwise, It is Ig- 
nored. It can be applied either before or after a 
contrast expansion, and Cdn be reapplied several 
times In succession. it produces a very clean 
looking Image, with little apparent ill effect on 
the boundaries. However, the cleaning process may 
possibly eliminate texture, which may he useful to 
the experienced eye In detecting anamolles. 

Expansions of limited ranges of gray values 
are moat effective. The histogram of the raw Image 
Is useful In defining the ranges to be expanded. 
Several successive expansions revealed all of the 
known anamolles In the Images under Investigation. 
The histogram of figure 3 is shown In figure 4, 
and two expansions are shown In figures 5 and 6. 

In figure 3, the values between 127 and 137 
have been expanded to full scale. Note that this 
makes the lesions labeled "e" and "f" In figure 3 
visible. Similarly, In figure b, the values be- 
acon 140 and 1)0 have been expanded to full scale, 
making the lesion labeled 'V visible. 


Figure 4, - Histogram of figure 3 
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Cone 1 us Iona 

A systematic Investigation of aevaral Image 
enhancement functlona with varying Input paraaetera 
applied to three typea of Images produced tha follow- 
ing conclusions: 


1. The application of certain digital Image en- 
hancement techniques aubatant lally Increaaea the 
visibility of anaaol lea In toaograaui of human 
anatomy. 

2. Image sharpening la of little value, apparently 
because of the low signal to nolae ratio and the 
small variations In the regions of intereat. 


Figure 5. - Contrast expansion of figure ). Pointers 
Indicate locations of lesions. 


3. Most of the smoothing techniques "smear out" too 
much fine structure. Selective averaging of aln- 
gular picture eleaenta produces a clearer Image, but 
may poaalbly obacurs texture which la useful In 
Identifying anaanlies. 

4. Contrast expansion la of greateat value. How- 
ever, its use must be guided by use of the gray- 
scale histogram. 
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Figure 6. - Contrast expansion of figure 3. Pointer 
Indicates location of lesion. 
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